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Abstract 

It is possible that there exist some interactions between dark energy (DE) and dark matter 
(DM), and a suitable interaction can alleviate the coincidence problem. Several phenomenological 
interacting forms are proposed and are fitted with observations in the literature. In this paper 
we investigate the possible interaction in a way independent of specific interacting forms by use 
of observational data (SNe, BAO, CMB and Hubble parameter). We divide the whole range of 
redshift into a few bins and set the interacting term 5{z) to be a constant in each redshift bin. 
We consider four parameterizations of the equation of state Wde for DE and find that 5{z) is likely 
to cross the non-interacting {5 = 0) and have an oscillation form. It suggests that to study the 
interaction between DE and DM, more general phenomenological forms of the interacting term 
should be considered. 
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I. INTRODUCTION 



It has been suggested from astronomical observations that the main components of our 
universe are dark matter (DM) and dark energy (DE). DM behaves hke the usual baryon 
matter and can form clusters, while DE is uniformly distributed in the whole universe, and 
it derives the universe to accelerating expand. Very ironically, we have known a little on 
DM and DE so far. The questions of what particles DM's are and what nature DE is remain 
open. From astronomical observations, however, some properties of DE can be deduced. 
For example, usually one characterizes DE with its equation of state Wde-, the ratio of the 
pressure to the energy density of DE; Wde is found very close to —1 from the observations. 
Therefore a natural candidate of DE is the well-known cosmological constant introduced by 
Einstein in 1917, for which the equation of state is exact —1. Although the cosmological 
constant is a beautiful and economic candidate, it suffers from some theoretical puzzles to be 
explained as currently observed DE. The theoretical difficulties (puzzles) are so-called fine- 
tuning problem and coincident problem (i.e., why energy densities of DE and DM happen 
to be of the same order today?) 

To avoid these problems, some dynamical DE models have also been proposed in the 
literature. The simplest dynamical DE model is a time-dependent scalar field. Based on 
different forms of the Lagrangian of scalar field, the scalar field models could be classified into 
quintessence, K-essence, phantom and quintom models. Furthermore, due to the ignorance 
for DM and DE, one is not sure whether there exists any direct interaction between DM and 
DE, at least no known symmetries prevent such interaction. Indeed, possible interactions 
between DM and DE have been intensively investigated in recent years. It has been shown 



that a suitable interaction can help to alleviate the coincidence problem 



l|-y] . Various 



interacting models have been studied 5l4l5l|. Several phenomenological interacting forms 



have been proposed and have been fitted with observations [16N22 |. Some recent discussions 
seemingly imply that the decaying of DM into DE is favored 23|, \2S[ by observations, which 
can make the coincidence problem more severe. However, most of those studies depend 
on the interacting forms, that is to say, those results are obtained by taking some special 
interacting terms. In other words, those studies are model dependent. Moreover, most of 
the models exclude the possibility of an oscillation interaction. And if the interaction exists, 
by fitting, one could only conclude that either DM decays to DE, or DE decays to DM. 
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In this paper, we are going to study the interaction in a way independent of the interacting 
form by observational data. To do that, we divide the whole redshift rarige into a few bins 



and the interacting term 6{z) is set to be a constant in each bin [25|, |27|. Clearly such 
study depends on DE models and the number of bins. We will study 3-6 bins cases with a 
constant Wde and try to get some common features of S{z). To see effect for different DE 
models, we will adopt four different parameterizations of Wdej^ith a preferable division of 
jins. We will fit the interacting models with the Union Snia [28], BAO [29], 9 Hubble data 
30l ] and the shift parameter R from WMAP5 31i- We obtain the best-fitted parameters and 
likelihoods by using the MCMC method. We find that 6{z) is likely to be oscillating and 
to cross the non-interacting {6 = 0) line. We also compare behaviors of r = p^/ p^e in the 
best-fitted models with those of corresponding models without interaction. In three cases 
of four parameterizations of Wde, the coincidence problem is alleviated, though DM decays 
into DE in some regions of redshift. 



II. METHODOLOGY 



We consider interacting models in a fiat FRW universe 



SH"^ = + Pb + Pde + Pd-n 



(1) 



where p^ and pf, are energy densities of radiation and baryon, respectively, and pde and pdm 
are energy densities of DE and DM, respectively. We have set the speed of light c = 1 and 
SttG = 1. The continuity equations for energy densities of the interacting DM and DE are 



Pdm + ^Hpdm = 
Pde + SH{1 + Wde) Pde = -SH6. 



(2) 



In some phenomenological models of interactioii, the interacting term 5 is always assumed to 
36 a function of pdm and pde, such as 5 = Xpdm [20, lllf , 5 = Xpde [24 132| or 5 = X{pdm + Pde) 
2|], thus the constraints resulting from observations will depend on the form of S. By fitting, 
if A = 0, it indicates that there does not exist interaction between DM and DE; if A > 0, it 
stands for the decay direction from DE to DM; while from DM to DE, if A < 0. However, 
obviously the way loses the possibility that 6 has an oscillating behavior. 
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To investigate such a possibility, in this paper we divide the whole redshift into four bins 
and set 5 to be a piecewise constant in each redshift bin 

5{Zn-l < Z< Zn) = 5n, > 1) (3) 

In our main analysis, we will set zq = 0, Zi = 0.2, Z2 = 0.5, z^ = 1.8 and Z4 = 1090. Also we 
will consider possible effect of the number of bins on the fitting results. 
From Eq. ([2]) we have 

Pdm{z) = pL(1 + - (1 + Zf [ -^^^dx, 



1+X 



p..W=p»,FW + fw£|^^|l^dx, (4) 

where F{z) = exp[JJ' ^^^i^^dx] and superscript represents the present value. As pb = 
p°(l + z)^, we can write pdm and pb together as pm- With the piecewise constant 6{z) we 
have an analytical form for pm 

Pm{Zn-l <Z<Zn) = [p'L - ^1 + J^^^' ~ ^^+^'>^^ + + + (5) 



J) — J) i J) 



where = + p^^. 

For energy density pde of DE, we will employ four different parameterizations of w^e as 
follows. 

I. Wde = —1- In that case, pde can be written analytically as 

Pde{Zn-l < Z < Zn) = Pde + ^ ^(^i ~ ^^(^ + + ^^^i ln(l + z) (6) 

i=l 

II. Wde = Wq. In this case, we have 

Pde{Zn-l < Z<Zn) = [pde + , , ^ 

I + Wq 



III. Wde = Wo + Wiz/{1 + z). In this case, F{z) in Eq. (j4]) has the form [33|, |34| 



F{z) = (1 + ;2)3(i+-o+«.i) exp(-^!^) (8) 

1 H" 

IV. Wde = '^0 + Wiz/{1 + zf' , In this case F{z) can be expressed as 

F(z) = (l + .)3(i+-)exp(^^^) (9) 
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For parameterizations III and IV, it is hard to get analytic forms of pde as Eqs. ([6]) and 
([7]). Now the Friedmann equation of the interacting models can be written as: 

E\z) = = + + pjml + p,jml (10) 



We now fit these four models with observations. The observational data to be used are 

the Baryon Acoustic Oscillation (BAO) data from SDSS DR7 



;he 307 Union SNIa data 



29|. the shift parameter R from WMAP5 3l|], and 9 data of the Hubble parameter B.(^z^ 



30l |. We obtain the best-fitted parameters by minimizing 

xL = Xln + xL ^X\^x\^{h- 0.742)70.0362 (11) 

where h = i^o/lOOkm ■ s'^ ■ Mpc" . We acquire the constraints by using the MCMC method. 
For 307 Union SNIa data, xL is defined as 

Xsn 2^ ^2(^A I ) 

where fith{z) = 51og]^o[(l + ^) Jo d-^/^i^)] + f^o, and po = 42.384 — Slog^^o^ is a nuisance 
parameter. One can expand Eq. (fT2|) as 



xln = A + 2poB + ^ilC 



where 



, _ [fJ'thiZi; /ip = 0) - flob{Zi)f 

We adopt the minimization of xL with respect to /io to replace xL 

Xsn Xsn,min JD / \-y 



In fact, it is equivalent to performing an uniform marginalization over /xq 36|. 
For the BAO data, one has 

2 [Dy(0.35)/Dv.(0.2) - 1.736]^ 

Xbao Q_Qg52 ^ ) 

where Dv{z) = [z/ H{z){^ dx/ H{x)fYl'^ . 



For the shift parameter, we take 

0.0192 



, _ 1.71)2 

Xr - n nin2 l"^^^ 



where i? = Jq' and 2, = 1090. 

And for the Hubble evolution data, we have 



^.^^j^mz^y^H^ (16) 

i=l 



Note that we also have used a Gaussian prior h = 0.742 ± 0.036 [37]. 



III. RESULTS 



Now we fit our models with the observations. As the data only give very weak constraint 
for z > 1.8, we fix 5(1.8 < z < 1090) = in our main analysis. To obtain the constraint 
for a specified parameter, we marginalize over all other parameters by using the MCMC 
method. In addition, in all computations, we demand that pde and pm keep positive in the 
range of z & (0, 1090). We do not decorrelate the constraints in the different redshift bins. 



The constraints of 6i are correlated. But in this way it ensures that the constraints o 



stained 



28 



35|. As 



for a given bin are confined to the exact redshift range of the bin, as discussed in 
the fitting results might depend on the divided method of redshift bins and the dark energy 
models, we will study two situations: 

A. different numbers of bins with a constant Wde', 

B. four different parameterization of Wde with a preferable division of bins. 



TABLE L The best-fitted values of 5i in z S (0, 1.8) for 3-6 bins. The values in () is the upper 
boundary of the redshift bin. 



Models 


Si 


82 63 64 


S5 




3 bins 


-0.32 (0.6) 


0.12 (1.2) 5.80 (1.8) 






4 bins 


-0.17 (0.45) 


-0.69 (0.9) 5.13 (1.35) -2.54 (1.8) 






5 bins 


0.05 (0.36) 


-1.09 (0.72) 0.16 (1.08) 10.12 (1.44) 


-10.32 (1.8) 




6 bins 


0.09 (0.3) 


-1.84 (0.6) 0.29 (0.9) 6.11 (1.2) 


1.08 (1.5) 


62.20 (1.8) 
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A. Effects of the number of bins 



At first, to see effects of the number and locations of the redshift bins, we divide the 
region of z E (0, 1.8) equally into 3, 4, 5 and 6 bins respectively and assume Wde to be a 
constant. The best-fitted 6i from the observations are shown in Table HI In the most cases 
the result prefers that 6{z) crosses the non-interaction line {6 = 0) around z = 0.5 and there 
is likely an oscillation of the interaction term 6{z). The corresponding errors of 6{z) are 
shown in the Fig. [H In 95% c.l., all results for 3-6 bins show that 6 is consistent with 
from the observations. But in 68% c.l., 6{z) is minus at the first bin {z ^ 0.5) in the 3 and 
4 bins cases and at the second bin (around z = 0.5) in the 5 and 6 bins cases . If one uses 
other three parameterizations of DE that introduced in the section II, one would get similar 
results, i.e., there is always a bin in which 6{z) departs from in 68% c.l. . As a result we 
may conclude from Fig. [1] that 
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FIG. 1: Constraints of 6{z) at 68% and 95% c.l. for 3, 4, 5 and 6 bins in z G (0, 1.8). The bins are 
equally divided. Tlie equation of state of DE is assumed as a constant. 5{z) is in tfie unit of 3H^. 
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1. With more bins, more finer structure of S can be resolved, e.g., for more than 3 bins 
oscillation behaviors of 6 appear. But for more bins the constraints of 6{z) in each bin from 
the observations will be weaker. 

2. The errors for z > 1 are much bigger than that for z G (0, 1). It is mainly due to the 
fact that there are much less data points in the large redshift region. 

It is also very likely that 6{z) crosses the 6 = line around z = 0.5. 



B. Effects of parameterizations of DE 



By considering these conclusions, in what follows we will divide the region of 2; G (0, 1.8) 
into three bins as: {zq = 0,zi = 0.2,^2 = 0.5,^3 = 1.8) [27|, which is the case adopted by 
most discussions in the literature and from which fine constraints of 6{z) could be obtained 
indeed. Four parameterizations of DE introduced in section II will be used. The best-fitted 
parameters and the constraints at 68% and 95% c.l. are shown in Table HTl The best-fitted 
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FIG. 2: The constraints of 5i at 68% and 95% c.l. . 5i is in the unit of 3/?q. 
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FIG. 3: The behaviors of r = pm/ Pde- Sohd curves are for the best-fitted models shown in Table 
im (with (54 = 0), while the dashed curves for the best-fitted models with 5{z) = 0. 

parameters for models with ^4 unfixed are also shown in Table HIl There are almost no 
differences between the models with ^4 fixed and unfixed. The corresponding 68% and 95% 
constraints are shown in Fig. |2l and the Fig. [3] shows the behaviors of the ratio r = pm/ Pde 
in the best-fitted models with interaction, compared with the cases without interaction. 
I. Wde = -1 

From Table [TTl and Fig. |2l it is obvious that 62 is negative, while 61 and ^3 are positive, 
which implies that 6{z) could cross the 6 = line at a recent time and the decay direction 
between DM and DE could be variable. When 6{z) crosses the non-interacting line, the 
effective equation of state of DE will cross the cosmological constant {w = —1) line. Thus 
DE behaves as a quintom fluid. As shown from Table HTl and Fig. HI 62 departs from beyond 
68% c.l. But models (111) and (IV), which have more degrees of freedom, show that it is 
still consistent with 6 = everywhere in 68% c.l. . In addition, it can be seen from Fig. [3] 
that the interaction can help to alleviate the coincidence problem in this case. 
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II. Wde = Wo 

For this parametrization, the situation is similar to the case of Wde = — 1- The sign of 
5{z) can be varied in the different bins. The possibihty of ^2 < is larger than 68%, but 
less than 95%. The coincidence problem is also alleviated in the best-fitted model. 

III. Wde = Wo + WiZ/{l + Z) 

In this case, there is still a downward departure of 82 from 0, but now the constraint is 
consistent with 5{z) = in 68% c.l. It looks from Fig. 3 that the coincidence problem could 
not be alleviated in this case. Note that to avoid a serious degeneracy, we have assumed a 
prior Q^rn < 0.37 here. 

IV. Wde = Wo + WiZ/ (1 + zY 

In this case, 5(yz) = is consistent with the observations in 68% c.l. There is still a 
possibility of crossing the non-interacting line. The coincidence problem can be alleviated. 

As expected, the resulting constraints are effected by parameterizations of Wde and di- 
visions of bins. But for all cases we have considered here we see that the interacting term 
prefers to have a behavior crossing the non-interacting line, and there might exist an oscil- 
lation 5{z) in the most cases. 

IV. CONCLUSION 

We have investigated the constraints of the interaction between DE and DM from the 
observational data. To make the constraints independent of specific interacting forms, we 
divide the whole redshift into four bins. In each bin 5{z) is set to be a constant 5j. First 
we have estimated effects of numbers of bins by consider 3-6 bins with a constant Wde-i from 
which we get some common features and choose a preferable division of bins: {zq = 0, = 
0.2,^2 = 0.5,2:3 = 1.8). For models of DE, we have adopted four parameterizations of Wde- 
The resulting constraints of 5j depend on these parameterizations. But there are also some 
common features of the interaction. The results are summarized as follows. 

1. The observational data prefer that 5{z) crosses the 5 — Q hne and has an oscillation 
behavior at a recent time, which implies that the decay direction can be variable. It is 
similar to the case that the equation of state of DE is likely to cross the cosmological 
constant {w = —1) line. For many well studied phenomenological interacting forms, such 
as S — Xpdm and S = X{pdm + Pde), the sign of S{z) is unchangeable. Our results raise the 
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possibility that S{z) can have different signs at the different times. It imphes that more 
general phenomenological forms of the interaction should be considered, if the interaction 
indeed exists. 

2. The constraints given from observations show a departure of 6{z) from beyond 68% 
c.l. for the Wde = — 1 and Wde = wq parameterizations. But for other two parameterizations 
of DE, Wde — Wo + Wiz/{1 + z) and Wde — wq + Wiz/{1 + zY which have more degrees 
of freedom, the constraints are consistent with 5{z) — To confirm the existence of the 
interaction, more observations and theoretical studies are needed. 

3. The coincidence problem can be alleviated in the three cases of four models, compared 
to corresponding ones without the interaction. The ratio r = p„J pde will evolve more rapidly 
(slowly) when d{z) < { S{z) > ) than the cases without the interaction. The decay of 
DE to DM {S{z) > 0) can alleviate the coincidence problem, while the decay of DM to DE 
{S{z) < 0) will make it more severe. Though S{z) is negative somewhere in the best- fitted 
models, its effects can be offset by that of the S{z) > regions and the period of r ~ can 
be longer than that of the corresponding non-interacting models. This way the coincidence 
problem is alleviated. 

Due to the ignorance on the properties of DE and DM, to study the dark sector inter- 
actions one always needs to assume some models of DE and DM. DM is always assumed 
as pressureless fluid, while there exist plenty variants of DE models. In the sense of phe- 
nomenology, the main difference among those models of DE is just the number of parameters. 
In our discussions, we have adopted four widely used parameterizations for the equation of 
state of DE in the literature. For different numbers of bins, there are also common results. 
Therefore our results are of some universality in some sense that 6{z) prefers to cross the 
(5 = line and have an oscillation behavior. In particular, this results indicate that if there 
does not exist any interaction between DM and DE, the model of DE should be paid spe- 
cial attention with an oscillating equation of state, because the oscillating behavior of the 
interacting form is mathematically equivalent to the case without interaction, but with an 
oscillating equation of state of DE. 
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EoS of DE 




h 




Wo 


Wl 


Si 


S2 


Ss 


^4 


I.w = -1 {64 = 0) 


ML 

best-fitted 


n 70+0.02+0.05 

'-'•'"5-0.02-0.06 

0.734 


n OS+0.04+U.08 
^•^°-0.06-0.12 

0.244 


{-1} 

{-1} 


{0} 
{0} 


n 1 0+0.22+0.42 

^•■'■'5-0.19-0.41 

0.127 


n 7Q+0.43+0.92 
^•'^-0.40-0.61 

-0.804 


n c«+0.45+0.86 
^•'■'"-0.59-0.90 

0.770 


{0} 

{0} 


64 unfixed 


best-fitted 


0.735 


0.246 


{-1} 


{0} 


0.133 


-0.810 


0.778 


-0.097 


TT / r rv\ 

11. W = Wq {04 = U) 


ML 
best-fitted 


„ -O+0.03+0.06 
^-0.02-0.05 

0.735 


„ „„_(_0.04+0.08 
U-^°-0.05-0.13 

0.204 


n 0/^+0.11+0.18 
U.OU_o.28-0.72 

-0.792 


{0} 
{0} 


n no+0. 28+0.47 

0.484 


„ otr+0.75+1.33 
— U.OO_o.47_o.77 

-1.545 


„ r-n+0.45+1.00 

U.OU_o.60-0.88 
1.107 


{0} 
{0} 


S4 unfixed 


best-fitted 


0.735 


0.218 


-0.797 


{0} 


0.518 


-1.558 


1.128 


-0.603 


III. {64 = 0) 
w = wq + wiz/ (1 -|- z) 


ML 
best-fitted 


n 70+0.03+0.05 
'J- '•-'-0.02-0.05 

0.733 


n 00+0.05+0.08 
^■^^-0.06-0.11 

0.242 


A Q7+0.22+0.38 
^•y' -0.29-0.75 

-0.632 


1 99+0.37+0.75 
-■■•^^-1.81-2.75 

-1.507 


n. 9Q+0.27+0.44 
^•^y-0.45-0.79 

0.574 


n r;Q+0-67+1.26 
^•'3^-0.64-1.04 

-1.802 


n t;i +0.59+1.15 
^•^■■■-0.60-1.03 

1.134 


{0} 
{0} 


64 unfixed 


best-fitted 


0.734 


0.232 


-0.627 


-1.533 


0.551 


-1.801 


1.131 


-0.388 


IV. (^4 = 0) 

W = Wq + Wiz/ (1 + z)^ 


ML 
best-fitted 


n 7Q+0.02+0.05 
^■'^-0.02-0.06 

0.733 


n Qn+0.04+0.08 

U-^^-O.05-0.14 

0.221 


n QQ+O.65+1.67 

-0.594 


9 Q+2.0+4.1 
^■^-12.2-38.3 

-2.115 


n 9Q+0.57+0.89 
^■^^-0.34-0.55 

0.524 


n no+O-52+0.86 
^•'-'^-0.82-1. 53 

-1.807 


n co+0.38+0.81 
^•^°-0.53-1.39 

1.288 


{0} 
{0} 


84 unfixed 


best-fitted 


0.733 


0.228 


-0.598 


-2.098 


0.540 


-1.805 


1.287 


-0.277 



